Electronic and transport properties of kinked graphene 
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Local curvature, or bending, of a graphene sheet is known to increase the chemical reactivity pre- 
senting an opportunity for templated chemical functionalization. Using first principles calculations 
based on density functional theory (DFT) we investigate the reduction of the reaction barrier for 
adsorption of atomic hydrogen at linear bends in graphene. We find a significant lowering (~ 15%) 
for realistic radii of curvature (~20 A), and that adsorption along the linear bend leads to a stable 
linear kink. We compute the electronic transport properties of individual and multiple kink-lines, 
and demonstrate how these act as efficient barriers for electron transport. In particular, two parallel 
kink-lines form a graphene pseudo-nanoribbon structure with a semi-metallic/semi-conducting elec- 
tronic structure closely related to the corresponding isolated ribbons; the ribbon band gap translates 
into a transport gap for transport across the kink lines. We finally consider pseudo-ribbon based het- 
erostructures, and propose that such structures present a novel approach for band gap engineering 
in nanostructured graphene. 



I. INTRODUCTION 

Nanostructures based on graphene have an enormous 
potential for applications. Especially in future elec- 
tronic devices compatible with and extending silicon- 
technology, due to the outstanding electronic transport 
properties of graphene^. However, it is crucial to mod- 
ify the semi-metallic electronic structure of graphene to 
exploit its full potential for many electronic applica- 
tions: a band gap can be introduced by nanostructuring 
graphene. 

A common approach towards engineering the elec- 
tronic structure is to form quasi-lD graphene in the 
form of nanoribbons (GNR/2-1 The electronic structure 
of GNRs depends on width, direction and edge struc- 
ture - all parameters which to some degree can be con- 
trolled. GNRs can be formed by etching^, by unzipping 
carbon nano-tubes (CNTs^P, or, ultimately be grown 
with atomic-scale precision using self-assembly of precur- 
sor molecules on metal substrates^. However, for elec- 
tronic applications this approach requires a structure- 
preserving means of releasing and transferring the struc- 
tures to an insulating substrate. Bonding of H or other 
species to graphene with large coverage opens an in- 
solating band gap at the adsorption sites due to sp 3 
hybridizatiorP. Periodically ordered clusters of adsorbed 
hydrogen can be formed on graphene in patterns dictated 
by the Moire lattice mismatch between graphene and the 
metal substrate, which opens a semi- conducting band 
gap^. Finally, regular perforations, known as a graphene 
anti-dot lattice (GAL ffi or a nanoscale mesh of holeS^El 
can have neck widths^^ down to 5 nm corresponding 
to band-gaps of the order of 1 eV^l 

Graphene consists entirely of surface atoms and is thus 
exceedingly sensitive to the surroundings. In particular, 
the van der Waals (vdW) interaction with the substrate is 
of importance. The substrate interactions, which make 
graphene cling to small features, may be exploited by 
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Figure 1. (Color online) (a) Smooth ripple-like structure 
where first and last six rows of carbon-dimers are surface- 
clamped regions with a separation of L= 90 A. Atomic hydro- 
gen is adsorbed at positions I- VIII. (b) The resulting kinked 
graphene structure after hydrogen is adsorbed in lines at the 
most reactive position (II) corresponding to the smallest local 
radius of curvature. The four kinks divide the structure into 
five sections, S1S5. 

manufacturing nanostructures in the substrate. Periodic 
steps in a Cu substrate has been used to induce "wrin- 
kles" or ripples in graphene with period and height on 
the order of 10 mxP^. Recently, Hicks et al^ demon- 
strated how arrays of ID large-bandgap, semiconduct- 
ing graphene nanoribbons corresponding to a width of 
~1.4 nm can be formed in graphene on a step-patterned 
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SiC substrate. The substrate interactions can clamp a 
graphene sheet while partly suspended across small holes, 
so that a pressure difference on the in- and outside lead 
to the formation of bubbles or "blisters" in the sheetP^. 
Also, linear folds, where the graphene sheet is bulging 
up from the substrate, have been induced for graphene 
suspended over trenches using heat treatment 1 - 1 . Thus, 
the sheet can obtain significant bends at certain places 
induced by the substrate interaction, substrate nanos- 
tructuring, and subsequent treatments^. Calculations 
by Low et alW^ showed how a sharp step of height 1 nm 
in a SiC substrate, comparable to experimental values^, 
can induce a linear bend in the graphene sheet with a 
radius of curvature down to around 1 nm. 

The ability to accurately control sharp local curvatures 
of graphene presents opportunities for strain-assisted 
modification of the local electronic structure and chem- 
ical reactivity in the graphene sheet, and may open a 
route to band gap engineered devicefl^HsHII] Very re- 
cently Wu et a/pD showed how graphene on a Si substrate 
decorated with Si02 nanoparticles induced local regions 
of strain and increased reactivity in a selective manner. 
Atomic hydrogen or other chemical species do not easily 
react with flat graphene when dosed from a single side^l 
However, at positions where there is a substantial local 
bending, rippling or strain o f the graphene sheet the re- 
activity changes significantly 5 23 . So far there has been 
only a few theoretical studies of the atomic geometry of 
hydrogenated ripple structures in unsupported, strain- 
induced, graphene rippleiP^H. However, to the best of 
our knowledge no studies have addressed the reactivity of 
bends or the transport through hydrogenated ripples, or 
discussed the possibility of stabilizing non-planar struc- 
tures by hydrogenation. 

In this paper we consider the reactivity of linear bends 
in a graphene sheet, and the electronic transport proper- 
ties of kinks resulting from hydrogenation of bends. Our 
starting point is the generic graphene structure shown in 
Fig. [l^i which is inspired by the experimental observa- 
tion of trench formation^. The bulging of this structure 
results from shortening the distance between two sepa- 
rated, clamped regions in the sheet. The remaining sec- 
tions of the paper are organized as follows: In Section [H] 
we describe our computational method and setup. In 
Section we present our results: First we describe the 
adsorption barriers for reacting with single atomic H on 
the graphene bend at positions with different local cur- 
vature (positions I- VIII in Fig. [lj,) . Then we show how 
a linear bend transforms into a kink when decorated by 
H along the most reactive (most curved) line (Fig. [T]d) , 
and present the electronic transmission through a single 
kink in Section IIII Bl The kink acts as an effective bar- 
rier with transmission depending on kink-angle, ip. In 
Section [ill C| we study how two parallel kinks lead to the 
formation of a pseudo-ribbon type electronic structure. 
Finally, in Section |III D| we demonstrate the opening of 
a transport gap for multi-kink systems, such as the one 
shown in Fig. [T]d, and Section IV gives our conclusions. 



II. SYSTEM SETUP 



The bend we consider in Fig. [IJi is created along the 
armchair direction by fixing the first and last six rows 
of carbon atoms and shortening the separation L, while 
the rest of the atoms are allowed to relax. A separa- 
tion of L = 90 A is chosen in order to obtain realistic 
curvatures^. We have first assessed the reactivity of the 
structure at positions with different local curvature, see 
positions I-VIII in Fig. [T^,. Subsequently we have re- 
laxed the structure where lines of hydrogen (Fig.[lJ)) have 
been placed at the points of lowest radius of curvature, 
i. e. the points of highest local reactivity. This particular 
system is meant to illustrate the potential of the hydro- 
gen adsorption mechanism, and to gain insight into the 
modification of the electronic properties due to the hy- 
drogen lines. In a corresponding experimental setup we 
can imagine placing graphene across a trench, which al- 
lows hydrogen adsorption on either side of the sheet. 



The atomic and electronic structure calculations 
are based on density functional theory (DFT) using 
the SIESTA^ code, and the PBE-GGA exchange- 
correlation 38 functional. We employ periodic boundary 
conditions (PBC) in the direction along the bend with 
a cell-width of four carbon dimers, and 10 Monkhorst- 
Pack fc-points. We use a mesh cut-off of 500 Ry through- 
out. When calculating reactivity in the form of reduced 
reaction barriers the unit-cell is chosen so that the dis- 
tance between single hydrogen atoms is larger than 8.5 A. 
This ensures low hydrogen-hydrogen interaction, which is 
known to impact reaction barriers^. In the total energy 
calculations of relaxed atomic geometries and reaction 
barriers we also use PBC transverse to the bend (5 fc- 
points). We use a TZP basis-set for hydrogen and a SZ 
basis for carbon, except in the reaction barrier calcula- 
tions where we compare calculations using a DZP and 
SZP basis for the four carbon atoms nearest to the hy- 
drogen. In the barrier determinations we furthermore 
use spin-polarized calculations, because of unpaired elec- 
trons. For the relaxed geometries a force tolerance of 
0.01 eV/A is used, and the final energies are corrected 
for basis set superposition error (BSSE). 



Based on the computed atomic and electronic struc- 
tures we subsequently use the TranSIESTA 40 method to 
calculate the electronic conductance per transverse-bend 
unit-cell width. To this end we attach semi-infinite flat 
graphene electrodes to each side of the selected kinks, e.g. 
replacing sections S± and S2 in Fig. [TJa by semi-infinite 
electrodes in order to calculate the transmission through 
the single kink separating Si and S2 ■ In the conductance 
calculations we employ a dense transverse fc-point grid of 
400 points. 
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Figure 2. (Color online) Calculated reaction barriers for hy- 
drogenation of bend graphene as a function of the local radius 
of curvature (II- VIII in Fig. [l^). Flat graphene (position I) 
has an infinite radius of curvature and is used to normalize the 
barriers. Calculations are spin-polarised and allow for atomic 
relaxation. 



III. RESULTS 
A. Adsorption barrier 

Adsorption of hydrogen on graphene involves a reac- 
tion barrier which needs to be overcome before the single 
hydrogen atom sticks to the graphene sheet. Several in- 
vestigations based on DFT calculations show that atomic 
hydrogen adsorbs on-top on flat graphene with a barrier 
abo ut 0.2 eV and binding energy in the range of 0.7- 
1 e\^2H3Ml Th us a m inimum kinetic energy for the first 
hydrogen to reactP^U is required in an out-of-equilibrium 
situation such as in an atomic bearrPH. Casolo et alW^ 
calculated the reaction barrier and adsorption energy for 
multiple hydrogen atoms on flat graphene. In agreement 
with other studies they found decreased barriers to stick- 
ing for the second H atom, compared to the barrier for 
adsorbing a single H atom on a clean surface^. 

Here, we first focus on the trends in the change in ad- 
sorption barrier as a function of the local curvature of the 
graphene sheet. To this end we have considered atomic 
hydrogen absorption at the on-top carbon positions at 
points with different curvature on the bent structure, see 
Fig. [iji (positions I- VIII) . The barrier is determined by 
calculating the total energy for each position of hydro- 
gen above graphene as the hydrogen is moved succes- 
sively closer to the graphene. Following the adsorption 
investigations on flat graphene by Ivanovskaya et alv^* we 
perform a relaxation in each step of the hydrogen-bonded 
carbon atom and its three nearest neighbours. Using the 
method described above we obtain a reaction barrier of 
0.22 eV on locally flat graphene. This is comparable to 
results obtained by several groups using DZP or plane 



wave basis sets and the PW91 

find that a SZP basis set for the relaxed carbon atoms 
yields a reduced barrier height of 0.18 eV (both basis 
sets with orbital range corresponding to an energy-shift 
of 0.01 eV). Hence, we use the SZP basis in the following 
reaction barrier calculations in order to lower computa- 
tion time. 

For the positions (I- VIII) we obtain the reaction bar- 
rier for adsorption of hydrogen as a function of the lo- 
cal radius of curvature (RoC) shown in Fig. [2j The 
second least curved position (VIII), resulting in a large 
RoC, reduces the barrier by roughly 3% compared to flat 
graphene (position I). The most curved position (II) in 
our considered structure has a minimum RoC of - 20 A 
resulting in a barrier reduction of roughly 16%. For com- 
parison, this RoC corresponds roughly to the radius of a 
(25,25) nanotube. Experiments show that Cgo molecules 
are more reactive than CNTs, which in turn are more 
reactive than graphite^. Thus a significant increase in 
reactivity is expected for the graphene with the linear 
bend, and a simple Arrhenius estimate using our data 
yields a factor of 3-4 at room temperature (300 K). We 
have also performed calculations using the less rigorous 
DFTB method 34 and obtained results in agreement with 
the trend in reaction barrier reduction obtained above. 

We may understand the reaction barrier and its change 
with curvature by considering the changes in carbon bond 
lengths. The barrier is due to the fact that the react- 
ing carbon atom has to be pulled out of the graphene 
plane, stretching the strong carbon-carbon bonds, when 
reacting with the incoming hydrogen atom. When the 
graphene sheet is curved the carbon atom is already 
slightly out of the plane, and thus the energy required 
to pull the atom further out of plane is decreased com- 
pared to flat graphene. The ortho- and para-locations in 
the graphene hexagon have been shown to be the pre- 
ferred locati ons for hydrogen adsorption in studies of 
flat graphend 22 | 35 | 39 l With this in mind as well as the 
curvature-related reduction of reaction barriers we con- 
clude that the considered system allows the adsorption 
of hydrogen atoms in single lines along armchair-edges. 
The kink in the atomic structure due to the sp 3 -binding 
of a single H makes the graphene curve even more in its 
vicinity which, in turn, preferentially lowers the barrier 
for absorption of a H along the linear bend. This suggests 
a mechanism where the hydrogen adsorption is propagat- 
ing and leads to the decoration of the entire linear bend 
turning it into a kink-line. It may be viewed as analogous 
to crack formation mechanisms, where the breaking of a 
bond increases the stress on neighbouring bonds - only 
in this case, the graphene is hydrogenated rather than 
broken or destroyed. 



B. Single kink 

Next, we examine the energetic and transport prop- 
erties of kink-lines in the armchair-direction. We first 
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Figure 3. (Color online) Electronic transmission through a 
single kink normalized with the pristine graphene transmis- 
sion (To) as a function of the kink angle, ip, for electrons 
(E > 0) and holes (E < 0). The arrow indicates the normalised 
transmission at the equilibrium angle determined from the to- 
tal energy calculations shown in the inset. 



consider a single kink with angle ip, e.g. between sections 
Si and 5*2 in Fig. [TJd. The kink-angle ip is varied in the 
range 0° . . . 90° , while the three nearest unit cells on each 
side of the kink are allowed to relax. The total energy per 
H is shown as a function of ip in the inset of Fig. [3j show- 
ing a minimum energy for <p«50°. This angle roughly 
corresponds to the angle in a sp 3 -configuration where 
2(^9 = 109.5°. The adsorption of H causes local changes 
in the geometry, i.e. only the carbon atoms very close 
to the kink are moved, while the remaining structure re- 
mains unperturbed. For this reason, the adsorption of 
hydrogen atoms can be considered as a process which 
locally pins the bend. 

The electron transmission per unit-cell width is lin- 
ear in energy for pristine graphene in a energy range 
around the charge neutrality point (E = 0), e.g. TqCcE. 
We find similarly that the calculated kink-transmission 
curves also are linear, and therefore we express the re- 
sults for the transmissions in terms of the roughly en- 
ergy independent ratio T/Tq = const. The kink breaks 
the electron- hole symmetry and we fit E>0 and E<0 
separately, as shown in Fig. [3] Larger kink-angles result 
in an increase in the overall transmission, which may be 
contributed to a better 7r-orbital overlap across the kink. 
For the equilibrium angle, ip = 50° , the ratio T/Tq is close 
to 0.17 in both regions (indicated by the arrow in Fig. [3]), 
corresponding to a transmission reduction of 83%. Thus, 
we see that the hydrogen-induced kinks in graphene can 
be used to form effective electron barriers. We now turn 
to the effect of multiple barriers and periodic kink struc- 
tures in order to examine resonant tunnelling phenomena 
and band gap formation. 



C. Two kinks 

Band structure calculations show that periodic nano- 
rippling of the graphene is not sufficient to create a 
band gap^ due to the low scattering by the elastic 
deformation 17 . In contrast, periodic arrangements of 
adsorbed hydr ogen can indeed induce a semiconducting 
band gap 1 24 * 26 1 . The electronic band structures of hydro- 
gen lines on flat graph ene have been examined by Cher- 
nozatonskii et al) 25 ^ 26 \ and recently also for nano-rippled 
graphene^. Here we show how two parallel kinks lead to 
a local electronic structure which resembles that of an iso- 
lated GNR between the kinks. Such structures could be 
produced experimentally using the techniques described 
by Pan et al^. Hydrogen terminated armchair GNRs 
are semiconducting but have a small energy gap when 
the width in atomic lines is N = 3L — 1, where L is an 
integer^. In Fig. [4] we compare the electronic bandstruc- 
ture for armchair GNR (aGNR) (left panel) to the elec- 
tronic transmission through two kinks separated by the 
corresponding aGNR width (right panel) . In the present 
case the initial width (or, kink separation) is N = 17 
atomic lines of carbon, which shows a semi-metallic be- 
haviour in the transmission spectrum with a small trans- 
port gap. In accordance with isolated aGNRs the next 
two widths N = 18, 19 are semi- conducting, while the 
last investigated width N = 20 is semi-metallic again (not 
shown) . The close correspondence between the electronic 
band structure for the GNR and the transmission gap for 
the double kink system allows us to consider the struc- 
ture between two kinks as a pseudo-ribbon. 

For the semi- conducting pseudo-ribbons transport 
gaps surrounded by van Hove-type ID behaviour are seen 
in the transmission functions (Fig. 4j right panel). The 
transport gaps, E gap — 0.4 eV and i? gap = 0.5 eV, are 
in reasonable agreement with power-law scaling of E g&v 
with width found for aGNRipS. We note that the pseudo- 
ribbon breaks electron-hole symmetry. For the TV = 18 
case a larger van Hove resonance is seen at the valence 
band-edge, while for N — 19 a larger resonance is seen 
at the conduction band-edge. There are small transmis- 
sion values within the electronic band gap due to leakage 
through the barriers, which we expect to introduce shifts 
in the energies between the real and pseudo aGNR. 



D. Multiple kinks 

In order to illustrate the behaviour of systems with 
more kinks we consider a system consisting of four 
hydrogen-induced kinks, as illustrated in Figs.fljs and [5] 
The sections S1/S5 are now replaced by left/right infinite 
lead structures, and the "top" S3 pseudo-ribbon is con- 
nected to the leads via the "side" S%, S4 pseudo-ribbons. 
We keep £2, S3 identical for simplicity, and determine 
the transmission across the kinks (in the z-direction in 
Fig. [5|, which is experimentally more feasible. We now 
investigate how the different sections influence the to- 
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Figure 4. (Color online) (Left) Band structures for H-passivated armchair ribbons with varying width, N. The ribbons are a 
zero-gap semi-conductor, a semi-conductor with band gap _Bn+i=0.4 eV, and a semi-conductor with band gap _En+2 = 0.5 eV 
for widths N, N + 1, N + 2, respectively. All widths are based on the number of carbon atom lines N = 17. The band gaps are 
indicated by arrows and highlighted. (Right) The electronic transmission functions for the corresponding pseudo-ribbons, i.e. 
across two parallel kinks of varied separation as shown in inset. 
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Figure 5. (Color online) Projected band structure and transmission through structures with multiple kinks. The top (section 
53) and connecting pseudo-ribbons (sections S2/S4) are varied in width, changing their electronic properties. The projected 
band structure along the pseudo-ribbons (i.e. k x ) is shown using filled circles with radii proportional to the weight on section 
S3. Semi-metallic (M) pseudo-ribbons corresponding to TV = 17 are shown in blue, while semi-conducting (SC) ribbons of 
width N + 1 are shown in red. Bandgaps of the top sections -E ga p and connecting sections iSSJjf are highlighted. These are 
£ga d p = 0.39 eV (b), ElU = 0.42 eV (c), and S^g = 0.35 eV, E B S ^ = 0.27 eV (d). The transmission T e (E) (per simple 
transverse unitcell) is determined across all kinks in the z-direction, and transmission gaps comparable to the bandgaps are 
observed in subfigures b and d. 



tal transport for the four possible combinations of semi- 
conducting (SC) and semi-metallic (M), corresponding to 
the pseudo-ribbon widths N, N + 1 used in Fig. [4] 

In order to analyse the transmission we single out the 
band structure projected on to the top section, 6*3 (ex- 
cluding carbon and hydrogen atoms at the kink), in the 
band structure along the pseudo-ribbon direction. The 



weight on S3 is represented by a circle of radius R n k, 

Rnk oc \^{n,k x ,i)\ 2 . (1) 

ies 3 

Here "J/ is the wave function at k x (k z — 0) with n (i) 
being the band (orbital) index. The obtained projected 
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band structures are shown in the left parts of each sub- 
figure in Fig. [5] Generally, some bands have no weight 
(no circles), while others have significant weight indicat- 
ing that there is little mixing between the orbitals from 
section S3 and other sections. 

In Figs. [5^ and we consider pseudo-ribbons 
with semi- metallic top regions, namely S2/S3/S4 being 
M/M/M, and SC/M/SC, respectively. For the all-metal 
pseudo-ribbons, M/M/M, an almost energy-independent 
transmission function is seen with a transmission close 
to that of the metallic double-kink in Section fill CI The 
SC/M/SC structure shows a transport-gap similar to 
that of the single SC pseudo-ribbon with van Hove res- 
onances, while the S'3-projectcd band structure reveals 
isolated metallic states within the gap. We note that the 
transmission at the resonances for the SC/M/SC struc- 
ture is larger than the corresponding M/M/M transmis- 
sion. For the case of semi-conducting top pseudo-ribbons 
in Fig. ^ and Fig. [5ji, we note that M/SC/M show a 
greatly reduced transport gap compared to the single 
pseudo-ribbon case (also, note the scaling of the trans- 
mission axis), while the SC/SC/SC structure shows a 
complete extinction of the transmission in the electronic 
gap, as expected. Generally, we find that the main be- 
haviour of the transmission is controlled by the connect- 
ing sections S^j&j, i.e. there is a good correspondence 
between the side section bandgaps and the trans- 

mission gaps. 

IV. CONCLUSIONS 

The presented investigations show that linear kink-line 
structures may form in graphene by reacting with atomic 
hydrogen along a linear bend in the sheet. The adsorp- 
tion barrier is lowered in the close vicinity of the bend 
which can be exploited to form the kink. In particular, we 
have shown that a radius of curvature of ~ 20 A reduces 
the hydrogen adsorption barrier by roughly 16% com- 
pared to H adsorption on pristine graphene. The calcu- 



lations suggest that once a single hydrogen atom has been 
adsorbed, the induced local kink and resulting increase in 
local curvature makes it easier for the following H to ad- 
sorb, thus creating a propagating kink formation. A full 
line of hydrogen atoms pins the structure and divides the 
electronic systems into different regions. We have shown 
that the electronic transmission through a single kink is 
reduced by 83% compared to pristine graphene, meaning 
that the kink-line acts as an efficient barrier for electron 
motion. We have demonstrated how two close-by paral- 
lel kinks form a pseudo graphene nanoribbon with similar 
behaviour of the electronic structure as for isolated nano- 
ribbons. The transmission function displays transport 
gap features corresponding to the isolated nano-ribbon 
band gaps. 

The present work thus suggests that it may be feasi- 
ble to template functional electronic nanostructures us- 
ing the conformation of graphene, e.g. to the substrate, 
and that this in turn induce changes in local reactivity. 
Our work clearly calls for extensions in a number of direc- 
tions. First of all more calculations are needed in order to 
investigate the kink-line propagation reaction proposed 
by our results. To this end it is important to include 
a realistic description of the actual substrate. It is also 
interesting to consider other adsorbate species, possibly 
introducing doping of the pseudo-ribbons and electronic 
gating. Finally, decoration and pinning of the edges of 
other geometries such as "bubbles" or "blisters" is of 
interest, e.g. in order to produce GAL-like structures^ 
without perforating the graphene sheet. 
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